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Using a straightforward extension of the analysis of Lieb and Wu, we derive a simple analytic 
form for the ground state energy of a one-dimensional Hubbard ring in the atomic limit. This result 
is valid for an arbitrary number of lattice sites L and electrons N < L. Furthermore, our analysis, 
including an application of the theory of stochastic matrices, provides insight into the degeneracy 
and spin properties of the ground states in the atomic limit. We give numerical results which 
illustrate how the atomic limit is approached. 
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I. INTRODUCTION 

For nearly fifty years, the Hubbard modeli has been 
used to describe many-body effects in solids, captur- 
ing the dominant competition between the delocalizing 
effects of the kinetic energy (with strength described 
by a hopping energy t) and the localizing effects of 
the electron-electron repulsion (with strength described 
by an on-site Coulomb energy U). Despite its simple 
form, the Hubbard model has provided significant in- 
sight into many-body properties of solids such as metal- 
insulator transitions, high-temperature superconductiv- 
ity, and magnetic states^ largely because of the accessi- 
bility of its analytic and numerical solutions. 

The analytic understanding of the Hubbard model 
stems primarily from the seminal work of Lieb and Wu 3 -^ 
who derived a "Bcthe anzatz" method^ for determining 
eigenvalues and eigenfunctions of the single band, one- 
dimensional Hubbard model with L lattice sites and ./V 
electrons. They obtained an explicit expression for the 
ground state energy for a half-filled system in the ther- 
modynamic limit (./V = L — > oo). In the intervening 
years, many authors have extended the Lieb-Wu anal- 
ysis. For example, Essler, Korepin, and Schoutens&i 
have shown that for even L and periodic boundary condi- 
tions, it is possible to find operators to extend the Bethe 
ansatz solutions and to find a complete set of eigenstates 
for the model. Lieb and Wu^ and Goldbaum^ have pre- 
sented some existence proofs for Bethe ansatz solutions 
for ground states. In addition, many authors have ap- 
plied the Lieb-Wu equations to particular systemsi 9 ' 10 : 11 

In the present paper, we reanalyze the Lieb-Wu solu- 
tions to the single band, one-dimensional Hubbard model 
for the case that there are L lattice sites (L < oo) with 
periodic boundary conditions and N electrons (N < L) 
in the limit U/t = u — > oo. This represents the strongly 
correlated regime of the system for which we find that 
the exact ground state energy can be expressed by the 
simple analytic form 
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To the best of our knowledge, this exact result has not 
appeared in the literature as the u — > oo ground state 
energy of the one-dimensional Hubbard model with pe- 
riodic boundary conditions for general values of L and 
N < L. On the other hand, there are some closely re- 
lated results ; 12 : 13 ' 14 : 15 and it is consistent with the results 
derived for the thermodynamic limit i 16 : 17 ' 18 The simple 
analytic form {l} is helpful for analysis and numerical 
studies of the one-dimensional Hubbard model and its 
derivation provides insight into the nature of the eigen- 
states of the one-dimensional Hubbard model. 

In Sec. (|TTJ we present the derivation of the ground 
state energy expression ([1]). The analysis, which requires 
a simple extension of the Lieb-Wu approach, helps to il- 
luminate properties of the model. In Sec. we show 
several examples of systems with fixed L and N, illus- 
trating how total energy eigenvalues converge from their 
finite to infinite u values. In Sec. (|IV[) we discuss the 
consequent properties of the ground state wavefunctions 
for this system and summarize our results in the context 
of the literature. 



II. DERIVATION OF GROUND STATE 
ENERGY 

The Hamiltonian of the Hubbard model isi 

7~i = Whop(i) + 7~lint(U) 

= ~ l C i,<r C r,« + U Y1 n qA n q,l ' ( 2 ) 
(q,r) q 

where cj a ( c q,a) creates (annihilates) an electron with 
spin a in the Wannier state localized at lattice site q, 
and (q, r) denotes that a sum is over nearest neighbor 
sites only. For the case of a one-dimensional system with 
periodic boundary conditions, the site index q (or r) takes 
values 1 < q < L and indices q and q + L are equivalent. 

Lieb and Wu considered solutions to the one- 
dimensional Hubbard model for a specified electron spin 
distribution with iVj and iVj_ indicating the total number 
of up and down z-component spin orientations, respec- 
tively, where N = iVj + iVj_. Applying the Bcthe ansatz, 
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the total energy eigenvalues are written in a form identi- 
cal to that for independent electrons^ 
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S(JVi,JV T ) = -2t5]cos*y 
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The so-called charge momenta kj, however, are not the 
wavevectors one finds in the absence of electron-electron 
interactions. To find the charge momenta requires solv- 
ing the set of coupled nonlinear equations (the Lieb-Wu 
equations) 
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Lkj {N l , 7V T ) = 2-kIj + 2 E tan_1 
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In this formulation, Lieb and Wu assume N < L and 
JVj_ < Nf. The parameter Ij = Ij^N^) is an integer (half- 
odd-integer) if iV| is even (odd), where 1 < j < N. The 
parameter J a = J a (N-\) is an integer (half-odd- integer) 
if iV T = N - N l is odd (even), where 1 < a < iVj.. The 
A's are a set of ordered, unequal real numbers Ai < A2 < 
... < \N r Details of the derivation of these equations are 
given in Yang's examination of a one-dimensional system 
with delta function interaction^ as well as a more recent 
review of the Hubbard model by Lieb and Wu.— 

In the present paper, we focus on the limit u — > 00, 
henceforth referred to as the atomic limit. This limit, 
which represents a highly correlated system, simplifies 
the mathematical properties of the Hubbard model con- 
siderably. In particular, since we are considering the case 
where N < L, it is reasonable to assume that the charge 
momenta for the ground state are real so that terms of the 
form (sin kj)/u in Eqs. (H} and ([5]) vanish in the atomic 
limit. This assumption is not valid if N > L since the 
ground state energy is linear in u, making the fc^'s nec- 
essarily complex^ As a result, Eq. ([5]) simplifies to 



27V tan" 1 ) = 2tt J a + 2 ^ tan" 1 

13=1 
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. (6) 

By substituting this expression into Eq. ((JJ) we obtain an 
equation for the charge momenta at u = 00, 
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In order to analyze the ground state, one possibility^ 
is to choose Ij and J a to be consecutive integers (or half- 



odd-intcgers) centered around the origin. We will con- 
sider other possibilities later. With this choice, if N is 
even, 



E J = °; 



(8) 
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otherwise, if TV is odd, 



Nl 
13=1 



2 ' 



(9) 



It thus follows that the charge momenta have the follow- 
ing relationships 



k j (N l ,N ] )={ 
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N even 



N odd. 
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For a given system of L sites and N < L electrons, we 
now show that the lowest energy obtained by substituting 
the charge momenta given in Eq. (|10p into the energy 
expression given in Eq. ([3]) is not necessarily the energy 
of the ground state. To do so, we present some examples 
where this is not the case. 

First, we consider a configuration with all electrons 
having the same spin (iVj, iVf) = (0, N). The eigenstates 
of the Hubbard model in this configuration can be repre- 
sented in terms of an antisymmetrized product of inde- 
pendent electron states in a Bloch basis, where the single 
particle energies are given by 



e,- = — 2fcos 



2nj 



(11) 



In this expression, the j's are a set of any L consecutive 
integers. The minimum energy of this system is thus 
obtained by summing N terms of the Bloch contributions 
(fTTj) . resulting in the analytic form for £ , m ; n (0,iV): 
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N odd, 



(12) 

where E g is the energy given in Eq. (JTJ) . The form of Eq. 
(fT2"| shows an interesting simplification: if N = L — 1 is 
odd, then £' min (0,iV) = -2t. Interestingly, Trugmar3 
proved the general result that the lowest bound energy of 
the ground state of the one dimensional Hubbard model 
in the atomic limit is exactly —2t(L — N), where L — N is 
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FIG. 1: Plots of the exact ground state energy given in Eq. 
([1]) (dashed line) and the lowest energy obtained using the 
Bethe ansatz charge momenta Eq. (|13|l in Eq. ([3]) (solid line) 
for a system with N — L — 1 = 4n + 1, as a function of n. 
The circles denote integer values of the variable n. 



the number of holes. Therefore it follows that the ground 
state energy of a one-hole system with an odd number of 
electrons must be E g = —It. Furthermore, there must 
exist precisely one state with that energy with maximal 
total spin. 

We continue our examination of a one-hole system with 
the special case2i oi N = L — 1 = 4n + 1 and iV[ = 2n, 
where n is an integer. In this case, Eq. (jTUJ) becomes 



kj(2n,2n + l) 
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where j is an integer. We find the possible energies of 
the system by summing over j in the formula given in 
Eq. ([3]); the lowest energy obtained is plotted in Fig. Q] 
as a function n. Inspection of this plot reveals that using 
the charge momenta given in Eq. (Tl"3"|) to evaluate Eq. 
© results in energies that are above the ground state 
energy of — 2t, except in the one-electron case (n = 0) 
and the thermodynamic limit (n = oo). The failure of 
this example for systems of finite size, indicates that Eq. 
(fT0|) should be reexamined. 

To address this problem, we recall that the Hubbard 
Hamiltonian obeys the commutation relations 



[n, S 2 ] = [H, S z ] = [H, S±] = 0, 



(14) 



where S 2 and S z are the operators for total spin and z- 
componcnt of spin, and S+ and 5_ are the spin raising 
and lowering operators. Because of these relations, it 
is possible to find the energy eigenstates of the Hubbard 
Hamiltonian which are also eigenstates of S 2 and S z with 
quantum numbers S and M s , respectively. 

The use of the charge momenta relations in Eq. (fT0|) 
does not take this spin symmetry into account. In or- 
der to include this symmetry in the Bethe ansatz for- 
mulation, we consider the relationships between systems 
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FIG. 2: Plots of the exact ground state energies using Eq. ([T]) 
(circles) and energies obtained using the unmodified Bethe 
ansatz solutions in Eqs. (|3J and (|1U|I (squares) at u = oo 
for L — 10 (top) and L — 15 (bottom) as a function of the 
chain density N/L. The dotted line connecting the circles is 
obtained by treating iV as a continuous variable. 



having different spin configurations. We assume that the 
energy eigenfunctions tp(Ni, N^) are known 



HV(^i,^t)=^V'(^,^ t ) ! 



(15) 



and are simultaneously eigenfunctions of S 2 with total 
spin quantum number S. Applying the spin-raising op- 
erator S+ to Eq. (|15[) . one of two things occurs: either 
the eigenfunction ip(Ni, N-f) is annihilated if it has mini- 
mal total spin quantum number S = ^ (iVf — iVj), or we 
obtain the eigenvalue equation 

Hi>(N i - 1, JV t + 1) = Ei>(N i - 1, JV T + 1), (16) 



where E is the same energy that appears in Eq. (|15p . 
Therefore, the set of eigenvectors and eigenvalues that 
solves Eq. (|16[) is a subset of the set of eigenvectors and 
eigenvalues that solves Eq. (fl5|) . This process can be 
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repeated to span all of the possible related spin con- 
figurations corresponding to the same energy eigenvalue 
E(M, M 1 ) = E(Ni,Nf), where the possible values of 
(M, M') are determined from limits of the raising op- 
erations to be < M < Ni, < M' < N, and 
M + M' = N. 

In determining the ground state energy of our system, 
we are not actually constructing the eigenstate wavefunc- 
tions, but searching for the minimum energy represented 
by Eq. ((3]) through the use of the charge momenta in 
Eq. (fTQ|) . The result of the spin raising operator analy- 
sis, suggests that Eq. (fTO)) should be modified to take the 
form: 



r 2^ 
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k j (N i ,N r ) = { 



M 



(17) 



where < M < N l . 

With this modified result, we now minimize the sys- 
tem energy in Eq. ([3]). We write the charge momenta 
appearing in Eq. (|3|) in a general form consistent with 
Eq. ([HD 
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where the j's are positive integers and jo is a real num- 
ber. With this choice of charge momenta kj, Eq. ([3]) can 
be summed as a geometric series resulting in the energy 
equation 
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This expression has a minimum when 

N + 1 
.7o = ~ — ■ 



(19) 
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Now we must check whether this mathematical mini- 
mum is consistent with the modified Bethe ansatz charge 
momenta kj(N±, iV|) given by Eq. (|17[) . First consider 
the case of ./V even. For this case Eq. (|2"D|) is a half-odd- 
integer and Ij(M) = j + jo must be a half-odd-integer, 
which is realized when M is odd. Since for all even 
N > 0, there is at least one choice of odd M in the 
range < M < N/2, this case is consistent with the 
modified Bethe ansatz solution. Now consider the case 
of N odd. For this case, Eq. (f20|) must be an integer and 
Ij{M) + M/(2N) = j + jo must be an integer, which is 
realized when AI = 0. To summarize all of these possibil- 
ities we conclude that the charge momenta corresponding 



TABLE I: Examples of minimum energies obtained using Eqs. 
© and (0 and particular choices of the J a parameters. The 
ground state energies E g were found using Eq. ([T]) and were 
verified by exact diagonalization. 
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to the ground state can be chosen using consecutive in- 
tegers j centered at the origin of the form 
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Using these charge momenta, and minimizing Eq. ([3]), 
we obtain the ground state energy expression in Eq. (JTJ) . 
This concludes our derivation. 

Additional confirmation of Eq. ([T]) was obtained by 
direct diagonalization of the Hubbard Hamiltonian in the 
u — ► co limit. Furthermore, it is easy to show that Eq. ([1]) 
is consistent with previously found lower bounds to the 
ground state energy^ In Fig. [5] a comparison between 
the lowest energies found using Eqs. Q and (jTUJ) , and the 
energy found using Eq. |T]) is shown for L = 10 and L ~ 
15 at various lattice densities. On this scale, the error 
introduced by using Eq. Q10p is small but non-trivial. 

In addition to the analysis presented above, there are 
other possible ways to find the ground state energies from 
the Lieb-Wu equations. Instead of restricting the J Q 's 
to be consecutive integers or half-odd-integers centered 
around the origin as assumed in the above derivation of 
Eq. (fTJ), it is also possible to find energies (including the 
ground state energy) using Eqs. © and ([7]) but with 
non-consecutive J Q 's. In Table U we list several examples 
of this situation. In all cases we have investigated, the 
ground state energies determined in this way are consis- 
tent with Eq. |T]). One might hope that such an approach 
could generate excited states in addition to the ground 
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FIG. 3: Plot of the lowest energies as a function of u for (top) 
TV = L = 5 and (bottom) TV = L = 6. A thick (red) line 
represents a maximal spin S = TV/2 state; a dotted (green) 
curve represents a 5 = (TV — 2) /2 state; a dashed (blue) curve 
represents a S = (TV — 4)/2 state. In the bottom plot, a solid 
(black) curve represents a S — (TV — 6)/2 = state. Symbols 
denote a doubly degenerate energy. Arrows mark the energies 
obtained by Eq. (JXJ) . 



states of the system. Unfortunately, there is no guarantee 
that the Bethe ansatz solutions are complete; in fact for 
even L it has been shown^il that extensions are needed 
to obtain a complete set of solutions. 



III. NUMERICAL EXAMPLES 

In order to understand the results presented here, we 
consider some specific examples. Throughout this sec- 
tion, we assume, without loss of generality, that TVj — 
TV| = (+1) if TV is even (odd). For ground states in the 
atomic limit, the probability of two electrons occupying 
the same lattice site vanishes. Therefore, there are only 
two configurations for a lattice site: vacant (henceforth 



FIG. 4: Plot of the lowest energies as a function of u for (top) 
TV — L — 1 = 5 and (bottom) TV = L — 1 = 6, using the same 
line (color) style as Figure [3] The degeneracy of the ground 
state in the atomic limit is given by Eqs. (|28[1 (top) and (|33[) 
(bottom). 



referred to as a hole) or singly occupied. In general, for a 
system of n holes, the ground state is highly degenerate, 
with a degeneracy denoted by dW. 

We begin our examination with the simplest case (half- 
filling) where the number of holes is zero. For any half- 
filled system in the atomic limit, the ground states have 
energy E g = 0. The degeneracy of the ground states is 
the number of possible ways of filling TV = L holes with 
TV indistinguishable spin- 1/2 particles 
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J\r T !JV;! 



(22) 



In Figure [31 the energies of these d'- -' states, that be- 
come the ground states in the atomic limit, are plotted 
as a function of u for two small half- filled systems: TV = 5 
and TV = 6. Results were obtained by exact diagonaliza- 
tion. In either case, the degeneracy of the ground state 
in the atomic limit is given by Eq. (|2"2"|) : for TV = 5, 
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FIG. 5: Plot of the lowest energies as a function of u for (top) 
N = L — 2 = 5 and (bottom) N — L — 2 = 6, using the same 
line (color) style as Figure [3] 



d (a ~> = 10; for N = 6, d (0 > = 20. As predicted by Eq. 
(|12p . in both cases, since E g = 0, one of the ground 
states has maximal total spin. This state is easily distin- 
guished by having an energy that is independent of the 
coupling parameter u. 

Next, we introduce one hole and examine its effect. 
In Figure [4] the lowest energies of two small one-hole 
systems, N = L — 1 = 5 and N = 6, are plotted as a 
function of u. According to Eq. ([1]), the energy of the 
ground states in the atomic limit for any one hole sys- 
tem is Eg — —It. In order to determine the degeneracy, 
we use perturbation theory, taking 7Yi n t as the unper- 
turbed Hamiltonian. The unperturbed ground states are 
degenerate and have zero energy 
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This unperturbed degeneracy D is given by the number 
of possible ways of filling L holes with N = L — 1 indis- 



tinguishable spin-1/2 particles 

(N + iy. 



D 



(24) 



Treating the hopping term 7Yhop as the perturbation, 
the first order correction to the ground state energy is 
the lowest eigenvalue of the matrix 



w« = {k\n hop (t)\£). 



(25) 



In the atomic limit, higher order corrections to the 
ground state energy vanish; the lowest energy eigenstates 
of the perturbation matrix in Eq. (|25p are the exact 
ground states in the atomic limit. Furthermore, we ob- 
serve that this perturbation matrix can be put in block 
diagonal form, by appropriate ordering of the states from 
Eq. (|23p ■ To prove this assertion we examine more closely 
these states. 

In Figurc[H]we show several unperturbed ground states 
defined in Eq. ([23|) for the case N = L — 1 = 7 and 
Ni = 3. With periodic boundary conditions, the state 
depicted in (5') can be obtained from arrangement (5) 
by repeatedly moving electrons one at a time to near- 
est neighbor sites with the provision that all intermedi- 
ate states contain no doubly occupied lattice sites; the 
same cannot be said for states (1) through (4). There- 
fore, states (5) and (5') are said to be connected to each 
other and disconnected from states (1) through (4). In 
fact, states (1) through (5) represent the five distinct 
spin configurations for that system. To indicate that the 
unperturbed ground states can be sorted into different 
spin configurations, we write them as \a,i a ), where the 
integer a denotes to which spin configuration the state 
belongs, and i a is an arbitrary state label. Since states 
in different spin configurations arc disconnected, 



(a' ,i' a ,\Hhop\a,i a ) = 0, if a ^ a'. 



(26) 



Therefore, the first order perturbation matrix is block 
diagonal, with the number of blocks being the number of 
distinct spin configurations for that system. 

In general, the number of blocks is a complicated func- 
tion that depends upon our choice of TV and L. We treat 
the simplest case, when N is odd, first. Then, it is easy 
to show that the number of connected arrangements in 
each spin configuration is 



C = NL = N(N + 1); 



(27) 



this expression is valid provided that N a ^ and N^/Ni 
and its inverse are not integers for N a > 1. With the 



assumption that N-t 
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1, this expression is neces- 



sarily valid for an odd number of electrons. Therefore, 
the number of distinct spin configurations is 
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the last equality will be shown below. 
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Returning to Eq. 
tion as 



25l) , wc factor the first-order correc- 
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-2t{k\H hop (-£) 



(29) 



the lowest eigenvalue of W(i) corresponds to the maxi- 
mum eigenvalue of the dimensionless operator H.hop(— \)- 
By appropriately choosing the unperturbed basis defined 
in Eq. (|2"3"1) . it can be shown that if N is odd, then each 
row and column of Whop (— \ ) consists of only two nonzero 
elements, which have the value |. Therefore, 



(fc|W ho p (-5) \t) > 0, (30) 
and the sum of the elements in every row or column is 1 



|fc) 
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and 



5>|Whop (-1) \£) = 1. 



(31) 
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FIG. 6: Possible one-dimensional spin arrangements for N = 
L — 1 = 7 and JVj = 3. Arrangements (1) through (4) each 
represent a different spin configuration for this system; ar- 
rangements (5) and (5') represent the fifth and final spin con- 
figuration. Arrangements (1) through (4) are disconnected 
from one another and from arrangements (5) and (5'); ar- 
rangements (5) and (5') are connected to one another. 



Nonncgative matrices that satisfy either Eq. (|3Tj) or Eq. 
(|32|) are known as stochastic matrices^ and are well stud- 
ied. These matrices describe the transitions of a Markov 
chain; their elements are the transition probabilities that 
a system will jump from one state to another. Equations 
([31) and (|32|) state that the total probability of transition 
is unity. By the Perron-Frobenius theorem^, the max- 
imal eigenvalue of these matrices is always 1, with the 
corresponding unnormalizcd cigenstate being the unity 
vector, whose elements are all 1. 

Since each block is a stochastic matrix, we conclude 
that the ground state is er^-fold degenerate. Further- 
more, from Eq. (|2"9"|) . the energy of the ground state is 
—2t, in agreement with Eq. (Q]). For the case N = L — 1 = 
5 and 7V T = N l + 1 = 3, Eq. (J35J) predicts that in the 
atomic limit, dP' = 2, which is confirmed in Figure 2] 
For small values of u > 0, the ground state energy is 
two-fold degenerate. As we increase u, an energy cross- 
ing occurs in the vicinity u ps 100 and a different state 
becomes the ground state. As u — > oo, this state con- 
verges to the energy of the lowest-energy maximal spin 
state, which is — 2t. In the atomic limit, the higher en- 
ergy states plotted in Figure 2] have energy — 1.95630t; 
this result corresponds to the case of n = 1 in Figure 
[TJ which was obtained by using the charge momenta in 
Eq. (fP3"| . This was identified as the ground state in the 
atomic limit in Ref. l2ll . 

For an even number of electrons, the situation is 
slightly more complicated, since Eq. (|2"T| is not valid for 
every spin configuration. Particular spin configurations 
may have additional periodicities that decrease the num- 
ber of connected states in that configuration. However, 
with our assumption that TVj = N±, we find that if N a 
is prime, then there is only one spin configuration which 
does not obey Eq. (f2"T)) . This configuration contains all 



the states defined in Eq. (|2"3"|) where every pair of electron 
neighbors have opposite spin (a one-hole antiferromag- 
netic chain). Obviously, the number of connected states 
in this spin configuration is 2L. Therefore, the degener- 
acy is 



D-2L 
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For the case oi N = L — 1 = 6 and N a = 3 (shown in 
Figure 13), cfW = 4. 

The case of two holes in the atomic limit is much more 
complex, and the methods used for one hole are not ap- 
plicable. In Figure the lowest energies of two small 
two-hole systems, N = L — 2 = 5 and N = 6, arc plotted 
as a function of u. In general, the low-lying energies in 
Figure [5] show a striking similarity to the one-hole ener- 
gies plotted in Figure 2J Again, we find that in the case 
N = L — 2 = 5, there is an exchange of ground states in 
the vicinity u w 100. 



IV. DISCUSSION AND CONCLUSIONS 

The form of the ground state energy given by Eq. ((T|) 
and its derivation provides additional insight into the na- 
ture of the eigenstates of the one-dimensional Hubbard 
model in the atomic limit. The form of the energy Eq. 
([1]) shows that there is an electron-hole symmetry in the 
ground state energy such that the energy of a system 
with N electrons is identical to a system with L — N 
electrons, corresponding to L — TV and N holes, respec- 
tively. In the course of deriving Eq. ([T]), we used Eq. 
(fT2"|) to show that for an odd number of electrons N , or 
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at half-filling (N = L), the maximal total electron spin 
configuration has the ground state energy given in Eq. 
(Q]) for any system size L > N. The conclusion that 
for an odd number of electrons in the atomic limit, ex- 
actly one of the ground states (not counting the trivial 
2<S l max + 1 = N + 1 degeneracy) has maximal total spin 
has also been discussed by previous authors^. If N is 
even (except for the case of half-filling) , the minimum 
energy of a state with maximal total spin is above the 
ground state energy, as shown in Eq. (fl~2|) . 

Numerical results illustrate the asymptotic behavior 
of the energies of the ground state and low-lying excited 
states as u — > oo . In all the cases presented, the absolute 
error per particle of the u = oo energy given in Eq. ([1]) 
to the exact energy of the ground state at u = 100 is less 
than 0.03i per particle. 

Due to its rich structure and relative simplicity, there 
is an impressive literature devoted to solutions of the one- 
dimensional Hubbard model. To the best of our knowl- 
edge, Eq. |T]) and the derivation presented here has not 
appeared in the previous literature. On the other hand, 
there are closely related works. For example, Kotrla 13 
extended the approach of Caspcrs and Iskei2 to con- 
sider the u — > oo limit of the one dimensional Hubbard 
model from the view point of enumeration of all possi- 
ble single occupancy states of the system. The analysis 
of the minimum energy configuration results in an ex- 
pression that is equivalent to our Eq. (P) although the 
explicit analytic form is not given. Another related re- 
sult is by Schadschncider— , who augments the electron 
hopping term of the original Hubbard model in Eq. @ 
with a bond-charge interaction with strength parameter 
X. When X = t, the number of doubly occupied sites be- 



comes a conserved quantity. The energies of the modified 
Hamiltonian are determined. In the limit asm oo, the 
ground state energy of that model is the same as Eq. |T]). 
However, this is for X = t,u — > oo, whereas our result 
is for X = 0, u — > oo. Recently, Kumar— considered the 
fixed boundary solutions of a one-dimensional Hubbard 
system in the infinite u limit. For this case, the Lieb- 
Wu analysis is not applicable and the energy spectrum is 
quite different. Kumar was able to find the analogue of 
Eq. {I} for the fixed boundary case. 

More detailed analysis has been devoted to case where 
L is even which for periodic boundary conditions allows 
for bipartite symmetry^ Essler et al&l derive a method 
for finding all of the energy eigenstates by augmenting 
the Bethe ansatz using generators associated with the 
SO (4) symmetry of the system. Their results are pre- 
sumably consistent with those in this paper, though they 
do not explicitly evaluate their equations in the u — > oo 
limit. Licb and Wu^ and Goldbaum^ prove the existence 
of ground state solutions to the Bethe ansatz equations 
for the restricted case of even N = L and odd N a = N/2, 
and show that the ground state is non-degenerate. It 
should be noted that in this case, the ground state is 
unique only for finite values of u; in general, at u — oo 
the ground state is degenerate. Figure [3] illustrates this 
for the case N = L = 6: for u < oo the ground state is 
nondegenerate and has S = 0; at u = oo, d° = 20. 

In summary, we have derived an expression for the 
ground state energy of a Hubbard ring in the atomic limit 
for even and odd integer L. This expression agrees with 
exact diagonalization energies obtained for several small 
systems, and is consistent with limiting results reported 
in the literature . 16 ' 20 
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